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ABSTRACT 


In this paper, we apply two theoretical turbulence models, DIA and the recent GISS 
model, to study properties of a turbulent channel flow. Both models provide a turbulent 
kinetic energy spectral function E(k) as the solution of a non-linear equation: the two 
models employ the same source function but different closures. The source function is 
characterized by a rate n g (k) which is derived from the complex eigenvalues of the 
Orr-Sommerfeld (OS) equation in which the basic flow is taken to be of a Poiseuille type. 
The 0 S equation is solved for a variety of Reynolds numbers cones ponding to available 
experimental data. A physical argument is presented whereby the central line velocity 
characterizing the basic flow, Uq, is not to be identified with the Uq appearing in the 
experimental Reynolds number. A renormalization is suggested which has the effect of 
yielding growth rates of magnitude comparable to those calculated by Orszag and Patera 
based on their study of a secondary instability. From the practical point of view, this 
lenormalization frees us from having to solve the rather time consuming equations 
describing the secondary instability. This point is discussed further in XII. In the present 
treatment, the shear plays only the role of a source of energy to feed the turbulence and not 
the possible additional role of an interaction between the shear of the mean flow and the 
eddy vorticity that would give rise to resonance effects when the shear is equal to or larger 
than the eddy vorticities. The inclusion of this possible resonance phenomenon, which is 
not expected to affect the large eddy behavior and thus the bulk properties, is left for a 
future study. The theoretical results are compared with two types of experimental data: a) 
turbulence bulk properties, Table IV and b) properties that depend strongly on the 
structure of the turbulence spectrum at low wave numbers (i.e., large eddies), Tables V and 
VI. The latter data are taken from recent experiments measuring the changes in the 
propagation of an electromagnetic wave through a turbulent channel flow. The 
fluctuations in the refractive index of the turbulent medium are thought to be due to 
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pressure fluctuations whose spectral function II(k) is contributed mostly by the interaction 
between the mean flow and the turbulent velocity. The spectrum n(k) must be computed 
as a function of the wave number k, the position in the channel x 2 , and the width of the 
channel A. The only existing analytical expression for n(k), due to Kraichnan, cannot be 
used in the present case because it applies to the case x 2 =0 and A=co, which corresponds to 
the case of a flat plate, not a finite channel. A general expression for II(k,x,„A) is derived 
heie for the first time and employed to calculate the fraction of incoherent radiation 
scattered out of a coherent beam. In Section XI, we treat anisotropy and show how to 
extend the previous results to include an arbitrary degree of anisotropy in the sizes of the 
eddies. We show that the theoretical one-dimensional spectra yield a better fit to the data 
for a degree of anisotropy (o « 4) that is within the range of experimental values. We also 
extend the expression for FI(k,x 2 ;A) to II(k,x 2 ;A,o) and compute the pressure fluctuations 
for different values of a. Similarly, we evaluate the fraction of electromagnetic energy 
scattered by an anisotropic turbulent flow and find a good fit to the laboratory data for a 
value of a % 4 - 6. 

Theoretical problems however remain which will require further study: among them, 
lack of backscatter (i.e., the transfer of energy from large to small wavenumbers) in the 
GISS model, possible resonance effects between the shear and eddy vorticity, behavior of 
the one dimensional spectral function at low wavenumbers, and the role of the secondary 
instability. These topics are now under investigation. 
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I. INTRODUCTION 


Experimental data on turbulent channel flow 1-6 can be used to test the validity of 
theoretical descriptions of turbulence. The latter can be broadly divided into two 
categories: numerical simulation of the Navier-Stokes equations and theoretical closure 
models. 

While direct numerical simulations 7,8 have successfully reproduced several 
experimental data, they are limited to low Reynolds numbers since the number of grid 
points required increases according to the 9/4 power of the Reynolds number and the 
number of time steps needed for an accurate simulation increases according to the 3/4 
power, yielding a rate of increase of Reynolds number to the third power. For high 
Reynolds number flows, the required number of grid points rapidly outstrips presently 
available computational facilities. To treat high Reynolds number turbulent flows, large 
eddy simulations make use of (empirical) subgrid scale models 8 


Theoretical Closure models can be broadly divided into two categories: single— point 
and two-point closure models. The most well-known among the former is the one 
originally proposed by Hanjalic and Launder (HL),° which proved successful in describing 
several types of shear flows. The HL model provides three coupled differential equations 
for the Reynolds stress tensor r 12 , the energy dissipation rate c, and the turbulent kinetic 
energy K, defined as the integral over all wavenumbers k of the turbulent energy spectral 
function E(k). The latter is known to satisfy the following equation 10-1 " (Ref. 9, Eq. 4.4: 
Ref. 10, Eqs. 4.35 or 4.33; Ref. 12, Eq. 15.27), 


k k k 00 

|r E ( k ) dk + s [ E 12 (k) dk + [ k 2 E(k) dk = -[ T(k) dk , 

J o J o J 0 J k 
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where s is tne shear, E I2 is the spectral energy function of the Reynolds stress tensor, and 
T(k) is the non-linear transfer among eddies of different wavenumbers. Eq. (1) can easily 
be rewritten to exhibit the energy dissipation rate, c, by extending the k integration to 
infinity and then subtracting the result from (1). 

The HL model proposes an equation for this quantity which must be solved 
simultaneously with (1) and with the equation for e. As previously stated, the HL model 
has proven very successful in the detailed description of several types of shear flows, its 
main drawback being the presence of six parameters that the model cannot determine. 

The HL model had actually been preceded by the work of Tchen 11 who suggested a 
physically interesting model based on the possibility of resonance interactions when the 
eddy vorticity is equal to or greater than the shear of the mean flow. The need for an 
equation for E 12 was bypassed, for the model suggested an expression for it in terms of E(k) 
and T(k) themselves. Although the model introduces two free parameters, it proved 
successful in predicting not only the existence of an inertial Kolomogoroff region but also 
the existence of a k ^ range that was actually verified experimentally. Tchen's model is 
physically attractive and it would be interesting to try and improve on it. 

In the early seventies. Leslie “ was the first to consider in detail the possible 
application of two-point closure models to turbulent channel flow. Upon realizing the 
unmanageable complexity of the DIA equations for the general case of shear flow,^ Leslie 
tried to develop a systematic program of simplifications for the two— point closure equations 
by restricting the analysis to the case of anisotropic but homogeneous flow (i.e., with 
constant shear) in the hope of deriving in a deductive, parameter free fashion, the empirical 
one— point closure relations of Hanjalic and Launder.^ By his own surmising, Leslie did not 
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succeed in his attempt. 


Since Leslie's work was not followed by other attempts to employ two-point closure 
models to turbulent channel flow, one is left today with either the successful one-point HL 
closure model which, however, contains six free parameters to be determined from 
experiments, with direct numerical simulations with their intrinsic Reynolds number 
limitations, or with a large eddy simulation which must rely on the use of a subgrid scale 
model. 

The model we are about to present contains no free parameters, employs (two) closure 
models and. within the goals it sets for itself, can be said to be relatively successful. The 
models used in this paper contain features that are common to the HL model, to Leslie's 
attempt, and to Tchen's physical approach In common with the latter is the fact that our 
work deals with only one equation, the one for the turbulent kinetic energy spectral 
function E(k). but at the same time it differs from Tchen's work in that E 12 (k) is not 
written phenomenologically in terms of E(k) and T(k). Rather, it is considered in the same 
spirit of the HL work, namely derived from an external relation, see IV. Our model has in 
common with Leslie's approach the treatment of the non-linear transfer terms, but it 
overcomes the difficulty that he had in treating anisotropy (caused by the mean flow 
stretching the eddies in the streamwise direction) by adopting a less formal but physically 
appealing approach originally proposed by Kraichnan 28,29 We first derive expressions for 
all of the quantities of interest, i.e., the one— dimensional energy spectrum, the amplitude of 
the pressure fluctuations, and the attenuation of a laser beam in the stretched physical 
system where the wavenumbers are denoted by k'. These functions depend on E(k ' ). the 
turbulent energy spectral function in the anisotropic system. However, since present 
turbulence models only allow us to calculate the energy spectrum, E(k), in isotropic 
systems, we have introduced an aspect ratio or anisotropy parameter, o. a measure of the 
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stretching of the eddies in the streamwise direction, such that k; = k,/a, k 2 - k 2 . and k 3 
= k 3 . By transforming from the k' to the k system, we then perform the integrations in k 
space. The final expressions will depend on the anisotropy parameter a which is not known 
a priori but must be determined through comparison with experimental data. By 
comparing predicted and experimental one-dimensional spectra, we found that the best fit 
was for 4 < a < S. The amplitude of the pressure fluctuations at the walls is most sensitive 
to the value of a (since the flow tends toward isotropy at midchannel) and a value of a in 
the range 4 to 6 was found to best fit the measured values. In the calculation of me 
attenuation of a laser beam propagating through a channel flow, we again found a best fit 
to the data for a in the range of 4 to 6. The consistency of these results lends credence to 
the methodology employed. 

As stated earlier, we shall adopt the DIA as well a 5 the recent G1SS model. 14 The 
latter has been recently tested against several types of turbulence, convection, grid 
turbulence, shear, etc. In all cases considered, the results were satisfactory. 

Notation: while in the DIA model the spectral energy function is denoted by E(k). in 
the GISS model use is made of the function F(k)= 2E(k). 

II. THE GISS MODEL 

The GISS model is based on a physical representation of the non-linear energy 
transfer in a way that while reminiscent of the work of Heisenberg, 1 '^ generalizes it in two 
important aspects, as described below. Assuming the turbulence is stationary, we begin by 
writing Eq. (1) as 


i 


where 


f(k) = 2 J |n g (k) + i/k 2 E(k) dk 


(3) 


is the rate of energy input into the wavenumber interval [0,k], which is partly dissipated by 
molecular viscosity and partly transferred by the non-linear term, T(k). The DIA and 
GISS models differ in the way they describe the function T(k). 

Equation (2) separates the action of the source, represented by c(k), from the action 
of viscosity and the non-linearity. Clearly, the separation is not meant to imply that dk) 
does not depend on the same ingredients that enter into the right hand side of the equation, 
it only helps in visualizing the physics of Eq. (1) more clearly. The problem then reduces 
to the evaluation of the quantity n g (k). the rate at which energy is being fed into the 
system. A discussion of how to compute this quantity will be given in Section IV. For the 
time being, we shall only note that the structure of n g (k) must be such that it depends on 

the source of energy for the turbulent flow, namely the shear itself. That this is indeed the 
case, will be shown in Eqs. (53)— (54). 


We begin by writing the transfer term T(k) as the product of a turbulent viscosity 
times a mean square vorticity, i.e., 


J T(k) dk = x/ t (k)y(k) . 

k 


y(k) = 


k 2 F(k) dk , 


* 


(4) 


so that the equation describing the energy balance in a turbulent flow becomes 


e(k) = [i/+ i/ t (k)] y(k). 


If one adopts the expressions 


^ t (^) = | dk , <(k) = f = constant , (6) 

k 

Eq. (5) reduces to the well known Heisenberg model that successfully reproduces the 
inertial subrange form of F(k) a k when t/ is greater than v. 

The QI$g model adopts the form of Eos. (4) and (M but neither of Ens. (6) . 1, 
extends both expressions in such a way that they reduce to (6) only in the inertial 
subrange. The generalization is carried out in two steps. First, the expression for ((k). 
rather than being taken to be a constant independent of k, is written as 


<(k) = F(k) (n g (k) + vk 2 ) dk , 


where i/k 2 cancels the negative viscosity effects included in n (k): the latter represents the 

rate at which energy is pumped into the system in the wavenumber interval between 0 and 

k. It is expected that for sufficiently large k, the integrand vanishes since n — -ek\ 

s 

Thereafter. t(k) remains constant. Then the Heisenberg approximation, t(k) = constant = 
f is valid. 
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The second generalization consists of writing the turbulent viscosity as 


"t ,k) = I dk ■ (S' 

k 

where n c (k) is a correlation time between eddies, to be determined by a closure. Equations 
(7) and (8) represent a formal generalization of the Heisenberg expressions. The essence of 
the GISS model is contained in a new closure, i.e. an expression relating n c (k) to F(k) and 
n s (k), the latter being considered given, either from an independent calculation (see Section 
IV) or self-consistentlv.^ 

Differentiating Eq. (5) with respect to k and making use of Eqs. (7) and (S), one 
obtains 


n s (k) + H^i =k2l/ t (k) - W 

The GISS model closure is represented by 

k 2 ^ t (k) = >n c (k) , (10) 

which, when inserted into Eq. (9), yields an expression for n (k) in terms of n (k) and y(k). 

c s * 

namely, 


27n c (k) = n s (k) + ( n 2(k) + 47y(k)]^. ( 11 ) 

Substituting (11) into (9) and using (8), one obtains a non-linear differential equation for 
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F(k). This is equivalent to solving a non-linear equation for the auxiliary function V(k) 
defined as V(k) = v(k) + 7 n 2 (k)/ 2 , i.e., 


dVjjil = 2ing(k) 

dk k 


(12a) 


37n c (k) = n g (k) + [n s 2 (k) + 6 7 V(k)]* . 


Since dy(k)/dk = k 2 F(k), the knowledge of V(k) yields the spectral function. F(k). The 
paramete, 7 is the only free parameter in the model and was found to be related to the 
Kolmogorotf constant Ko by the expression 


7 = 



(13) 


One can easily show that the Heisenberg model is indeed contained in the G1SS 
model. For eddies sufficiently removed (in k space) from the source, which acts primarily 
at low wavenumbers, the time scale characterizing the source, n' 1 . is no longer the 
dominant one. The eddy dynamics are governed primarily by the local break-up 
mechanism characterized by a time scale given by the vorticity, so that when y* >> n , n 
^ y , and Eq. (8) reduces to the Heisenberg form given by Eq. (6), at least for a power law 
F(k). At the same time, far from the source, the function f(k) saturates, i.e., <(k) = ( = 
constant. Eq. (6). This shows that the GISS model encompasses the inertial subrange. On 
the other hand, for eddies with sizes of the order of the dimensions of the system itself, one 
may expect that the dominant time scale would be closer to that of the source i.e. n >> 

j s 

' ' n c ** n s' cleany f(k) is not constant. Under these conditions, one cannot expect a 
universal form for the energy spectral function, since F(k) now depends on the specific form 
of the stirring mechanism. 
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In conclusion, while the Heisenberg mode! yields a universal F(k), the GISS model 
yields an F(k) which is a universal function of n s (k). 

III. THE DIA MODEL 

Perhaps the most well known of the presently available theories to describe fully 
developed turbulence is Kraichnan's direct interaction approximation 18 (DIA) which in 
turn has given rise to other theories which bear the same spirit. We shall not discuss here 
the eddy damped quasi— normal Markovian *,RDQNM) approximation, which is 
phenomenological in nature, since the eddy correlation time scale (which within the DIA is 
determined by solving the integral equation satisfied by the infinitesimal response function) 
must be chosen using external inputs. A very complete description of EDQNM, its 
successes and limitations in describing "universal" properties can be found in Ref. 17. 

Since the DIA has been described in detail elsewhere, 12 we shall present only the basic 

equations for the turbulent energy spectral function. Although the DIA is a well 

understood approximation to the non-linear transfer terms, 18 it can be regarded as a fully 

deterministic theory without free parameters. Since the DIA formalism has in the past 

been applied primarily to describe those properties of turbulence that do not depend on the 

specific nature of the source function, experimental data concerning bulk properties could 

not be dealt with. To include them, the DIA was recently applied with good results to a 

set of model equations with a source function appropriate to high Rayleigh number 
19 

convection. Perhaps the main drawback in the application of the DIA formalism to 
specific cases of interest has been the rather intimidating nature of the equations describing 
the turbulent energy spectral function, E(k). Moreover, with the presence of the 
infinitesimal response function G(k), one must in fact solve two coupled integral equations 
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for E(k) and G(k). These equations may be written 19 




A 

Q(k,t-s) = 2 7r | J kqp b(k,q,p) dq dp 
s 

x [ J ds' G(k,s-s ' )Q(q,t-s ' )Q(p,t-s ' ) 

-oo 


ds' G(q,t-s')Q(p,t-s')Q(k,s-s') 


(14) 


-oo 


and 


(9 1 f f r 

~ n S ( k )J G(k,t-s) = -2tt J j kqp b(k,q,p) dq dp J ds' G(q,t-s')Q(p,t-s')G(k,s-s') 

s 

+ ( 15 ) 


where Q(k,t-«) - < Uj(k,t)uj(-k,s) >, u(k,t) is the Fourier component of the turbulent 
velocity , the angular brackets denote a realization average, t and s are time variables, the 
energy spectrum is given by E(k)=4 jrk 2 Q(k,0), b(k,q,p) = (q/k)(xy+z 3 ) with x, y, and z 
the cosines of the angles opposite k, q, and p, respectively, 6(x) is the Dirac delta function, 
and the A over the wavenumber integrals indicates that the region of integration is 

restricted to a subdomain in which k, q, and p form a triangle. Once the growth rate n (k) 

s 

is specified, Eqs. (14) and (15) can be solved for the energy spectrum and the infinitesimal 
response function. 
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IV. THE RATE n s (k) 


Neither the DIA nor the GISS model can be expected to fix the functional form of 
n s (k) which must be provided by considerations other than the ones that led to the 
equations for the kinetic energy spectral function F(k). Prescribing an n (k) is physically 
equivalent to prescribing an equation or a formula for the energy rate c(k) since 

r k 

f (k)=2 J (n g (k) +j/k 2 ) E(k) dk , (2) 

00 00 

f = f(oo) = 2 f (n g (k)+i/k 2 ) E(k) dk = 2v f k 2 E(k) dk . 
o o 

Our model can thus be considered a two equation model, one for the kinetic energy K. 

00 

K = | E(k) dk , (16) 


and the second equation for t(k). For a successful application of our model, the 
identification of the physically correct function n g (k) is clearly of critical importance. To 
that end, let us first discuss its physical meaning. From Eq. (9), we derive that 

n g (k) = k 2 i/ t (k) - n d (k) 

k 

v k > = k!E < k > dk ■ o) 
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i.e., the rate n g is the difference between two rates: the rate k 2 i/ t (k) at which energy, under 
the action of the non-linear interactions, cascades to all wavenumbers larger than the one 
considered (in fact, ^(k) is an integral from k to infinity), minus the rate at which energy 
is injected at k from all wavenumbers less than the one considered (since the the vorticity 
y(k) is an integral from 0 to k). 

Since n g represents a difference of rates characterizing the fully developed turbulent 
flow, i.e., long after the transition from the laminar state has occurred, it cannot a priori be 
identified with the instability function n g _T that characterized the transition between 
laminarity (L) and turbulence (T), even though the processes they represent are physically 
equivalent. To derive n g T , one has a well defined mathematical formalism, the stability 
theory. To derive n g (k) or, equivalently, e(k), one would have to construct another 
equation thus making the model much more complex. We have already shown 14 that 
physical arguments can be very helpful. Some general considerations are in order. First, it 
is known experimentally that, for example, in thermal convection the large scale structures 
that one observes at the transition do persist in the turbulent phase, i.e., their structure 
survives the strongly diffusive and shearing action of a turbulent flow. This result is 
perhaps not unexpected since the large scale structures have the longest lifetimes of all the 
eddies and also because their structure is affected primarily by the source rather than by 
the non-linear transfer interactions. Stated differently, since the largest eddies cannot 
originate from even larger ones, their sole source of growth is the source itself. Second, 
from a mathematical point of view it would clearly be greatly advantageous if one could 
employ even if partially, the well established mathematical framework of stability theory 
to gain information about the form of the function n (k). Here we want to make a clear 
distinction between the shape of the function n g (k) and its amplitude. We shall propose, 
and try to justify, that the former can be arrived at by the use of the Orr-Sommerfeld 
equation, while the latter can be arrived at only by providing a way to account for the 
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presence of a turbulent flow. The latter is in fact likely to renormalize in a significant way 
the amplitudes of the parameters characterizing the laminar regime. This in turn 
translates into a renormalization of the O— S equation itself. The success of the calculations 
presented in this paper will be seen to depend to a large extent on the proper 
implementation of this renormalization procedure. Since the latter depends on the specific 
problem at hand, it cannot be formulated in a universal fashion. Each physical problem 
brings in its own characteristic features. From the point of view of trying to understand a 
complex phenomenon like turbulence, we consider this to be an advantage, for the method 
requires an understanding of the renormalization that the most prominent features of the 
laminar flow have undergone. The use of an O-S type of equation to determine the 
functional form of the rate n g is proposed here because such a function has features of 
almost universal character that are bound to be sufficiently well described by such an 
approach. First, consider the shape of n g (k). For wavenumbers less than k 0 = 1/L, where 
L is the geometrical dimension of the system under consideration, there cannot be am- 
ercing and so n g must be less than or equal to zero. On the other hand, for large values of 
k, i.e., when one deals with small eddies, the dominant mechanism is kinematic viscosity 
which contributes a factor — «/k 2 , i.e., the function n g must become negative at some large 
value of k. In Ref. 14 it was shown that the GISS model requires that for large k, n ~ -vk~ 
quite independently of its behavior at low wavenumber, i.e., of the specific mechanism that 
feeds energy into the turbulent regime. One may therefore conclude that the general shape 
of this function must be of the form shown in Fig. (1). In our experience with different 
types of turbulence, i.e., grid turbulence, thermal convection, and shear, we have indeed 
verined that the physical n g (k) has the form of Fig. (1). For example, in the case of grid 
turbulence where both the kinetic energy spectral function E(k) and the non-linear transfer 
term T(k) have been measured experimentally, one can derive n (k) directly from the data 

o 

since from Eqs. (2) and (3) it follows that 
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2 n s (k) = "£[£) • O') 

The functions E(k) and T(k) are presented in Figs. (14) and (16) of Ref. 14. As one can 
see, the ratio does indeed have the shape of Fig. (1). 

A second, parallel argument as to the validity of an O-S type equation as a guide to 
the functional form of n g can be seen by using a procedure First suggested by Synge 20 
whereby the O-S equation 20-23 is formally rewritten so as to exhibit the instability 
function or growth rate n g (we shall omit the superscripts L-T). The result is (see Eq. 53 
below) 

n g (k) = Ar 12 (k) - t/B(k) , (18) 

which shows, as expected on physical grounds, that n g is composed of two terms: a source 
term, proportional to the shear r 12 , and a sink proportional to the viscosity u. The O-S 
equation further predicts that the shear peaks near the walls (see Fig. 4.21 of Ref. 23), thus 
implying that in that region energy is extracted from the mean flow and fed into 
turbulence. Experimentally, it is known (see Fig. 5.5 of Ref. 24) that the main source of 
energy production, i.e., of c(k), occurs precisely in that region. Thus, the O-S equation 
predicts correctly the physically important feature of a region of instability. There, the 
main physical process is intrinsically the same as the one that characterized the transition 
from laminarity to turbulence. That feature has thus survived even in the presence of 
turbulence, which one may hope can be accounted for by a process of renormalization for 
which we shall propose two methods and show that they yield very similar results. 


\\ e choose a coordinate system in which the mean flow is in the x— direction with one 
wall of the channel at y = — D and the other at y = D. Consider the well-known laminar 
Poiseuille profile, 24 with y- y/D, 
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(19) 


U(y) 


Di 

2 v 


[t 2 ] (*-/> • 


where Ap /( is the pressure drop along the channel of width 2D (we have taken unit 
density). In terms of the bulk velocity U m defined as the integral across the channel of 
U(y) divided by the width of the channel, 2D, i.e., 


m ZuT 


U D 
d rn 

m ~ v 


Eq. (19) becomes, where L stands for linear, 


V(y) = Uj (1 -/) , Ui = |u m , 


L 3 


and Uq is the value of U(y) at midchannel, y- 0. 


( 20 ) 


( 21 ) 


As one can see, in the linear regime the pressure drop is proportional to U . This is 
no longer valid in the turbulent regime, where the relation is quadratic rather than linear. 
We shall propose to incorporate this physical property into (19) as a way to renormalize for 
the presence of turbulence. In order to do so, it is convenient to introduce the friction 
coefficient A defined as^ 


A = 


I u 2 

8 u m 



( 22 ) 


so that Eq. (19) becomes 
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U <>-> = re7 U m O'/)' 


(23) 


While this new form does not introduce any new physics at the laminar level, being merely 
a rewriting in terms of different variables, it does allow us to incorporate turbulence effects 
if an additional relation for the quantity A is introduced, for example the one proposed by 
Blasius, 25 i.e., 


3A = fC 1 / 4 . 
m 


(24) 


Eq. (23) can then be written as 


U(y) = U*(l-/), (25) 

R 

where U 0 is the renormalized midchannel velocity. Had we assumed that Eq. (21) were 
valid in the turbulent regime, we would have 


U 


L _ 3 V 

0 “ 2T) 



(26) 


which shows that the renormalization procedure from Eq. (21) to Eq. (25) is equivalent to 
taking 


rrL 


0 ~ u » = 51 R m u 0 


(27) 


It remains to relate R 

m 


to the experimental R 


exp’ 



( 28 ) 
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where L 0 is the measured center line (i.e., midchannel) velocity. The relationship between 
R m and R exn is give " b .' 4 ' 2 ° 


R exp - R m “ 2 53 R r ’ R r = ^ • 


where the ratio R ex p/R r can be expressed alternatively as 3 ’ 4 ’^ (*=0.41, B=5.1 


j: In R r + B, 9.26S R° y ° 89 , 25 + 3.5 * 10 3 R . 


A second approach for renormalizing the laminar profile can be devised by considering that 
the experimentally measured mean velocity field, l T (y), does indeed retain some vestiges of 
the laminar profile in the region near the walls where the energy is being produced. This 
can be seen by recalling the experimental form of U(y) as given for example in Ref. 25. i.e.. 


U(y)«U r R r ( 1 - 1 / 1 ) R r (l-|/|)<l0 (3la) 

U(y)=U r (|lnR r (l-|^|) + 5.lj R J O-|/l)>10, (31b) 

and Fig. (5.5) of Ref. 24 which shows that the maximum production occurs at the point 
where the linear regime changes to the logarithmic one. To renormalize the profile (19), we 
must eliminate the pressure gradient in terms of some other physical characteristic of the 
turbulent regime. We shall proceed as follows. Since both the profile (19) and the y 
experimental one given by (31) are linear in the region where energy is being produced 
most efficiently, we shall require that the slope of the two functions be the same as y goes 
to -D so as to assure that the two are identical in that region. This implies that Ap /( = 
U^/D , which in turn implies that 
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(32) 


«(y)-u5(i-/), v*-^R». 

In this case we note that the renormalization for Uq is given by 

R ■ 

,,L i; R v r tT L 

L 0^0 = 3n IT - U o • (33) 

m 

In Table I, we present some numerical results that show how the two renormalizations. Eq. 

(27) and (33) yield similar results (r,= R;J/ 4 /72 , r 2 = R2/3R m ). Considering the 
empirical nature of some of the relations employed in the previous derivation, the 
agreement of the results calculated using the two methods may be considered quite good. 

It is probably not a coincidence that for R ex p = 10 4 i '«e get a renormalization factor 
that is verj close to the one first obtained by Reynolds and Tiederman“^ more than twentv 

V J » 

years ago and recently reconfirmed.^ Both studies point out that the use of the 

experimental form of U ex p(y>R exp ) in the O-S equation yields stable solutions, i.e., n g <0. 

f/it is assumed that the Reynolds number appearing explicitly in the O-S equation, Rq g. 

which is in effect a viscosity, is to be the same as the R gxp entering U exp (y,R e ). On the 

other hand, if R Q _g is left as a free parameter, unstable solutions (n>0) can be obtained if 

s 

R o _s-( 10 - 1 5) R exp. (34) 

It is as if the effect of turbulence is equivalent to a decrease of the molecular viscosity. At 
first this may sound incorrect for we know that when dealing with the mean flow equations, 
the molecular viscosity is increased by the effect of turbulence and, by extrapolation, one 
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might expect that the same would be true when considering not the equation for the rr. ..w 
flow but for the turbulent flow itself. The only phenomenon that comes to n. t nd that could 
play the role of an internal source of energy, thus a lowering of the viscosity, is backseat ♦er. 
a mechanism whereby the energy in the high wavenumber regions does not cascade to till 
higher wavenumbers (in the direction of the increasing importance of viscosity) but s 
scattered back to lower wavenumbers, thus in effect playing against viscosity. 


In conclusion, the rate n g entering in (3) is not the rate that characterizes the 
transition from laminarity to turbulence. Rather, it represents the rate at which energy is 
pumped into turbulence in such a waj as to satisfy (9), which represents the balance 
between the rate at which such an energy is put into the system and the non-linear terms, 
k 2 " t > th at either transfer it to higher wavenumbers or feed the same wavenumber interval 
with energy from lower wavenumbers, yCkjn" 1 ^). The quantity n $ represents an 
instability, but one in the presence of turbulence. In the case of channel flow, a region can 
be isolated near the walls where there is such an instability in the sense that the energy is 
being extracted from the mean flow and given to the turbulence. In that region, the 
velocity profile is linear. We have tried to extract information about the rate at which 
such energy is transferred from one type of flow to the other by using the O-S equation in 
which, however, the strength of the mean flow was adjusted to match the experimental one. 
In that sense, the 0— S equation must be viewed more as a phenomenological tool rather 
than the exact equation that it really is when one deals with the very different problem of 
the transition between laminarity and turbulence. 

There are, however, two other ways of looking at the rate n . One could renounce anv 

b ‘ 

attempt at deriving n g from a mathematical equation, but rather consider it as a 
phenomenological function of the type 
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n g = n g (k,a,b) , 


(35) 

where the analytical form can be taken to match Fig. 1, and where the parameters a and b 
can be Fixed by demanding that the total kinetic energy K and the rate ( derived from the 
turbulence model match the experimental values. By adopting this procedure, one 
transforms the model into one that, rather than predicting the turbulent quantities 
becomes a tool to assess the correctness of the non-linear transfer terms, i.e., the validity of 
the adopted closure model. This method was successfully adopted in Ref. 14 in the case of 
grid turbulence, where the predicted transfer term T(k) was found to compare rather well 
with the measured value. (See Fig. 16 of Ref. 14.) 

Another possible way of using the present model of turbulence is by way of inverting 
the model and consider tne basic equations (9)— (11) not as a tool to derive the spectral 
function F(k), and thus the turbulent quantities, but rather as an equation for n g itself. In 
fact using the experimentally measured one-dimensional spectrum Ej(k), one can construct 
the three dimensional E(k), i.e., 


2E(k) = 



(3G) 


the turbulent viscosity v (k) 

V 


"t< k > 



(37) 


and Finally the rate 7n c (k)=k'V t (k). Inserting these relations into (9) yields for n tin- 
expression 
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(38) 


n S ( k ) = 


^ 2k 4 I,(k) -I 2 (k) 

W k=lf(k) 


where 


I,(k) 


Mki. 


1 - 


EOT 



f^dk, 


k 


(39a) 


I 2 (k) = k 3 Ej(k) 


k dEn 

eot dlT 





dE^ 



-r 15 



dk , 


k 


0 


(39b) 


which can be easily computed once E,(k) is known experimentally. This procedure yields 

valuable information on the rate at which energy is being drained from the mean flow into 

ttnbulence. At the same time, the k dependence of n^ would help in identifying, or at least 

narrowing down, the possible types of tnstaoiiUi.es that might be candidates. Because of 

the integration from k to infinity in ( 39<- } , it is clear that a successful retrieval of n from 

s 

the experimental data can lie carried out only when Ej(k) is known over a large 
wavenumber interval. 


V. SOLUTION OF THE Ok Ft ? 'MMERFELD EQUATION 


We now proceed to show now the ;;ov-.,l 
the well known Orr-Sommrvf ia equation/ 0 " 


i ratt n . car; be ompu.ed from the solution of 
-23 

w nose i.erwation is briefly sketched below. 


mostly to firm up the notation and units employed. 


We begin with the Navier-Stokes equations for the total velocity v and write the 
linearized equations for Uj, where as before v.= 11+ Uj. Separating the other variables into 
mean and fluctuating parts, the latter being indicated by a prime, we obtain with l'. = 

no 1 


du 

3T 




(40) 

(41) 

(42) 


where u, v. and w are the x, y, and z components, respectively, of the fluctuating velocity 
and where d/dt = d/dl + U d/dx. In deriving the above equations, the mean flow- 
equations in the linearized limit were also used to eliminate the average pressure. The 
average density is denoted by p. Finally, the incompressibility condition is given by 


du dv dw 


= 0 . 


(43) 


Differentiate Eq. (40) with respect to x and Eq. (41) with respect to z, add the resulting 
equations and use Eq. (43) to obtain an equation containing only the velocity v. 
Differentiate this result with respect to y. Next, operate with ^/dx 2 +^/az 2 on Eq. (42). 
Using the resulting two equations to eliminate the pressure term, we deiive ( " = d 2 /dy 2 ) 

^V 2 - U "(y)^j v(x,y,z) = i/V 4 v(x,y,z) . (44) 
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Now assume that v has the form 


v (x,y,z) = v(y) exp [i(k x x -h k,z - k x ct)] . 


(45) 


With the notation 


k 2 = k 2 + k 2 , 
j. x z ’ 


Eq. (44) becomes the well known Orr-Sommerfeld equation,^ 


(46) 


v 

nr 

X 


|l-kj] 2 v(y) = (U(vH)[^-kJ 


V (y)-U"(y)v(y) 


(47) 


It is convenient to write this equation in terms of dimensionless quantities. Measuring 
lengths in terms of the channel half— width, D, and velocities in terms of the laminar 
centerline mean velocity, Uq, we introduce the following dimensionless variables: 


a=k x D, a=k x D, y/D, V=\J /uj, *= c/uj . (48) 

Eq. (47) then becomes 


i^f^- W-'-W = [*> - c] < - a»)»(> - v "(^ , 


(49) 


where 3= djdy The mean flow at this point is still arbitrary. For a Poiseiulle 
flow, is used and the boundary conditions v = d\jdy— 0 are imposed at 

/= * 1- If further define R L by oR L = qR l , we obtain 
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(50) 


r l = 


Vx 


K 


k;) j 


= R. cos d > , 

L 


taen the solution of Eq. (49) yields for a given R L the complex eigenvalue c as a function of 
o. Actually, for a given a, there is a large number of discrete complex eigenvalues c 
However, only one of the c’s has a positive imaginary part for some range of a which leads 
to growth of the instability (see Eq. 45). Defining the growth rate n = k .Im(c) (where Im 
stands for the imaginary part of the enclosed function), the solution of Eq. (49) yields 


n g = J 5 ~ o Im( c) cos <t> vs. Or, ( 51 ) 

for a given value of R. One can easily see that the maximum value of n (the most 

s 

dangerous mode) corresponds to taking cos <j> = 1, i.e., k =0. To compare our results 

z 

with experimental data, we need to change the velocity normalization from the undisturbed 
centerline mean velocity, T, T {j, to the experimental centerline mean velocity, U 0 . This is 
done by multiplying Eq. (51) by U 0 /Uq = R /R . Hence, the desired growth rate is 

hJ vs. k ± D, where n* = ^ , ( 52 ) 


for different values of R or, equivalently, of R = ILD/k 

L exp u 

VI. DETERMINATION OF k 2 

y 

Equation (52) is not yet the final result, for it gives n as a function of k . However. 

s ± 

since homogeneity has been assumed in the application of the DIA and in the construction 
of the GISS model, they require the knowledge of the growth rate as a function of the total 
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wavenumber, k. Channel flow is inhomogeneous in the direction perpendicular to the 
channel walls (the y-direction) so that an application of the model equations requires the 
construction of a wavenumber in the y-direction. i.e., k = k (k ). The application of a 
model which assumes homogeneity to a manifestly inhomogeneous flow may be criticized 
on first principles, yet we have learned from Leslie’s work 12 that the problem is 
analytically a very difficult one and trade-offs must be made to simplify it in order to 
make it solvable. We have kept as much as possible of the essential physics in the model 
intact, consistent with the necessity of keeping the problem solvable. One of these 
trade-offs is our inability to describe the variation in the details of the flow in the 
\— direction. All of our calculated bulk properties represent an average of these properties 
over the inhomogeneous direction. We therefore proceed to find a suitable relationship 
between k y and k^. Following the original work of Synge, 20 we multiply the O-S equation 
(47) from the left by v and integrate the result over y from -D to +D. Separating the 
real and imaginary parts, one obtains for the growth rate, n g = k x Im(c), the expression 

klm(Q) I 2 + 2k 2 I 2 + k 4 I 2 
S I?+k 2 I l I? + k 2 I 2 ’ ^ 

whe r e 
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In deriving Eq. (54), we have extended Eq. (45) to u and w and have used Eq. (43) in the 
form dv/dy = -i(k u + k w). 

A Z 


Equations (53)— (55) have a simple physical interpretation. The first term is a 
"production term", since it provides the interaction between the shear s = U'(y) in the 
mean flow and the Reynolds stress tensor: it gives a positive contribution to n , see Eq. (1 ). 
The second term is proportional to the viscosity, v, and is always negative. It must clearly 
dominate at very large wavenumbers where the only remaining physical mechanism is 
kinematic viscosity. As we know, the latter enters the Navier— Stokes equations in the form 
-i/k i 2 . On that basis, we shall therefore identify the coefficient of v with k 2 , i.e., 




+ k 2 = 
± 


J2 

*2 


+ 2k 2 I 2 + k*I 2 


I? 


+ kjl? 


which, after some rearrangements, yields the desired result 


(56) 



*2 + 

+ kjl l ’ 


(57) 


which we shall use to generate the wavenumber k from the solutions of the O-S equation. 

J 

For the fastest growing mode (cos 0=1), (57) becomes 



i 

' * * 

(u 3) 2 u + v ^ 2 v) dy 

i 


i 

* * * 

( u u + v v ) d^ 

- i 


(58) 
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This result is intuitively appealing since k 2 is found to be the average of the operator 
d 2 /dy 2 over the flow in the y-direction. 

VII. THE RESULTS FOR n (k) AND k 

S v 7 y 

We have solved the O-S equation (49) for four values of R . : 1.23, 2.86, 3.08. and 

4 ex P 

(x 10 ). The first three values correspond to the experimental values of Laufer 1 and 
Hussain and Reynolds. ^ In Table II, we present the values of n g (k) in units of n*, Eq. (52). 
as a function of the dimensionless wavenumber kD and in Table III we give the values of k 
vs. k^. In Fig.(2), we plot n g /n* vs. kD for different values of the Reynolds number. 

VIII. SOLUTIONS OF THE TURBULENCE EQUATIONS 

The function n g (k) was used to solve both the DIA and the GISS models, Eqs. (12) 

and (14)— (15). In each case, we computed several quantities of interest that we discuss 
below. 

1) Turbulent energy spectral function. F(k). In Fig. (3) we plot 

F(k)/F* vs. kD , F*= U 2 D = (iP/D) R 2 xp , (59) 

calculated using the GISS model. In Fig. (4), we compare the GISS and DIA results. 

2) The energy f(k). In Fig. (5) we plot the quantity 

f(k)/ f * vs. kD , e* = Ujj/D = (i/ 3 /D 4 ) R| xp , (60) 

where f(k) is defined in (5) and (7). (The physical units of e are erg g - ^ sec - ^.) 
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o» 


3) The turbulent viscosity, y (k). In Fig. (6) we plot 


v fi) !»* y s- kD, v* = U 0 D = v R exp , (61) 

where ^(k) is the turbulent viscosity defined in (8). The largest value of the turbulent 
viscosity is attained at k=k 0 , where k 0 is the smallest allowed wavenumber. Taking the 
limit k-k 0 in Eq. (5), we obtain 


^(ko) = 


n s (k 0 ) 

l - 2 


(62) 


4) One dimensional spectra. Both Laufer 1 and Hussain and Reynolds^ present their 
experimental results in terms of the one dimensional spectral energy function defined in 
terms of the three dimensional F(k) as follows (see Ref. 10, E- s. 3.72, 3.4S, 3.47) 


E 



F(k) (l-k?/k 2 ) dk. 


(63) 


We have computed E,(kj) for the case corresponding to the experimental condition in the 
Hussain and Reynolds paper, i.e., R exp = 28600, U Q = 1350 cm sec -1 , D = 3.18 cm. The 
results are presented in Fig. (7). 


On the other hand, Laufer 


1 


presents his experimental data in terms of the function 


F (p( n ) = U^ E i( k i) [J^Ei(k,) dk, 


(64) 
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where 2jm = k,U 0 . The comparison between Laufer's measured spectrum (R = 30S00. 
u 0 = "28 cm sec , D = 6.35 cm ) and the theoretical results is presented in Fig. (8). 

5) Turbulent velocities and scales of turbulence. Laufer 1 and Hussain and Reynolds 4 
also provide experimental data for other quantities of interest, namely 


a) turbulent energy, 


< u' 



dk , 


(65) 


b) turbulent velocity, 


<uj>* = <uj>* = <uM 
A y z 


r o ^ 

<u-> 

3 


* 


( 66 ) 


c) Taylor microscale A , 


K 


= f E,(k,)dk, [J”k?E,(k,)dk, 


-1 


(67) 


d) turbulence macroscale, A x , 


A 



Ei 



E(k) dk 



e) Kolmogoroff scale, 1 0 , 


lo = (^ 3 A) 1/ ' • 


( 68 ) 


(69) 
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In Table IV, we present <u 2 > in units of Ujj, <u;>* in units of U 0 , u.=u (k ) in 

x u t V V' 

units of U 0 D, e in units of U$/D, and A x , A x , and 1 0 in units of D. Experimental data are 
included when available. 

IX. PRESSURE FLUCTUATIONS IN A TURBULENT CHANNEL FLOW 

We begin by considering that in an incompressible fluid, the pressure p(x) is given by 
the solution of the Poisson's equation 


V 2 p(x)=-„ 




ckj dxj 


(70) 


where v(x) is the total velocity, p is the constant density, and the summation convention 
has been employed. If, as usual, one splits v. into a fluctuating part Uj and a mean field U-. 
the right hand side of (70) is seen to contain three contributions: a term in UjUj, one in the 
form UjUj, and one in UjUj. Following Kraichnan, 28,29 we shall assume that the cross 
term is much smaller than the other two so that from the solution of (70), one can 
construct the fluctuating pressure, i.e., 


^r- = J <* 3 k [n T _ T (k> + n TH| (k)] $ J d 3 k n;k) , (~i) 

where T-T indicates the contribution arising from the turbulence— turbulence interaction 
and T-M indicates the one arising from the interaction between the turbulence and the 
mean flow. It is a known fact, which we have verified numerically, that 


n_ >> n 

T— M T— T 


( 72 ) 
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The expression for ri T _ T was first evaluated by Batchelor. 30 The result is 

4* n T _ T (k) = || E(k-)E(k-k') | k S i n ^f| 4 k'Mk'sinftW . (73) 

The expression for ri T _ M is much more difficult to compute. To treat the physical problem 
under consideration, we need an expression that depends not only on k but also on the 
coordinate variable across the channel of finite width A, i.e. 

^T_y(k) "* • (74) 

The only expression for I1 T _ M available in the literature is one due to Kraichnan 29 (his Eq. 
5.20), which yields the pressure at the lower boundary of a semi-infinite medium, (a flat 
plate) i.e., 


• (75) 

Since we plan to study the propagation of an electromagnetic beam in a channel flow, 
Kraichan's expression is not applicable to our case. We shall therefore derive the 
expression for II T _ M (k,x 2 ) = II^k^A). 

Since in the case of a channel flow, two directions, say x, and x 3 , in the plane of the 
mean flow can be considered to be homogeneous, one can perform a Fourier transform of 
(70) on the variables x 3 , and u. Eq. (70) then becomes 


where k 2 = k^ + k 2 , and 


P(* 2 iM = (2 t) 3 / 2 J p(x,t)e *( k i x i +k 3 x 3 aA )dx 1 dx 3 dt . 


Eq. (76) is a second order, inhomogeneous differential equation whose solution can be 
written as 


x 2 

2k P (x 2 ,k,w) = J dx£ e~ K(x 2~ x 2)T( x ^/c,a;) - 

a 


x 2 

[ dx 2 e K ' X2 ~ x ^T(x 2 ,K,u) . 


(7S) 


The two constants of integration a and b must be determined by imposing the boundary 
conditions 




x 2 =0,A 


= 0 , 


(79) 


since the fluid is bounded by two walls at x 2 = 0 and Xj = A = 2D. The final result 


is 


U x 2 

P( x 2 ’M = | d x 2 g 2 (x 2 ,x^;K)T(x^) dx 2 gl (x 2 ,x$;,c)T(x$) , (SO) 


where, for the sake of simplicity, we have explicitly written only the x 2 dependence of T. 
i.e., T(x 2 ) = T(x 2 ,k,o>), and where we have defined two function gj and g 2 as 


Kg,(x 2 ,x£;*) = sinh k(x 2 -x^) 


(81) 


/t ?2( x 2’ x 2'^) ~ s i nh «A^ cosk ^(^~ x 2 ) • 


(Si 
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By taking the limit A - « and x 2 - 0, (80) reduces to Eq. (3.9) of Kraichnan. 29 (In tlu* 
case x 2 =0 and nonzero A, Eq. (80) does not reduce to Eq. 3.11 of Kraichnan 29 since the 
latter is missing a factor of 2 in front of the second term.) Eq. (80) is the basic ingredient 
of our calculations since it gives the value of the fluctuating pressure at any given point in 
a channel of full width A=2D. From (80) we obtain, with a slight change of notation, 

A A 

| p(x 2 ,K,o>) | 2 = | dx | dy g 2 (x 2 ,x;K)g 2 (x 2 ,y;K)<T*(x, /t , (i ,)T(y.«.u;)> 

A x, 

r r* 

-2 dx dy gi(x 2 ,y;K)g 2 (x 2 ,x;«)<T*(x,K,tj)T(y,/s,u;)> 
o 

x 2 *2 

+ J d x j dy gi(x 2 ,x;K)g,(x 2 ,y;«)<T (x,K,u>)T(y, *,«;)> , (S3) 

o o 

where the star denotes complex conjugation and the angular brackets denote an ensemble 
average. 

As is well known, the largest contribution to T(x 2 ) comes from the interaction of the 
turbulent field with the mean flow, i.e. 

T(x,t) = 2 /is(x 2 )tj^-2 , (84) 

where the shear s(x 2 ) is given by 
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( 85 ) 


so that 


<T(x.t)T(x',t')> = S(x 2 ,x 2 ,x 1 -x / 1 ,x 3 -x$,t-t') 

4/) S(x 2 )s(x 2 ) f^22( x 2’ X 2’Ci)C3iO i (86) 

where it has been explicitly indicated that in the directions x, and x 3 the fluid is 
homogeneous, R^ is the two point velocity correlation function, = x; - x,, and 
C 3 = x 3 ~ x 3 . Taking now 

R 22 ( x 2 ,x 2 ,^ 1 ,^3,t) = (2 t) 3 / 2 J d 2 «du;e*^ R^x^x^k,^) (87) 

and using an analogous expansion for S, we have 

S(x 2 ,x 2 ,/c,u>) — 4 p" s(x 2 )s(x 2 ) Kj R^x^x^,*,^) . ( 88 ) 

The relation between <T*(x,*u;)T(y,*,a/)> in (83) and the quantity S(x 2 ,x 2 ,«,^) is 29 
<T (x,n,u)T(y,K,uj)> = ( 2 x) 3 / 2 S(x 2 ,x£,k,u>). Using the reality of R 22 (x,t), we write 

i 00 

R 22 (X 2 ,x;,K,aO = ]£ C0S k 2 n) ( x 2~ x 2) MM - (89) 

n=-oo 

wheie k^ ^ =2jn/A. Eqs. (86)-(89) are substituted in (83) and the result is integrated 
over all u\ Using the relations 
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(90) 


00 

J Rj 2 (M du = (2ir) 2 $ 22 (k) 
—00 



where E(k) is the velocity spectrum for homogeneous turbulence, and where 

<p 2 (x 2 )> = | | p(x 2 ,k,u>) | 2 d 2 /tdw , (92) 


we obtain, after a series of lengthy integrations over the variables x and y in (83), the final 
result 
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with tfx) the Dirac delta function. To be consistent with the spirit of the model adopted in 
this paper, the velocity profile used in performing the integrations that lead to ( 93 ) was 
l (x 2 ) = s 0 x 2 (l - x 2 / A), where s 0 =r 0 / v. This gives for the shear 

s(x 2 ) = s 0 [l-^x 2 j. ( 97 ) 

However, it is clear that the above formalism is valid for any shear. 

Using the spectral functions E(k) derived from the DIA and GISS models, we have 

computed the pressure spectral function n(k) and then integrated over all wavevectors k. 

The resulting pressures, in units of are presented in Figs. ( 9 ) and ( 10 ). As one can 

see, for R exp = 5000, the calculated value of the pressure at the wall is 6.31, while 
31 

experiments and direct numerical simulation 32 yield a lower value, i.e. 3 . 22 . the 
discrepancy being due most likely to anisotropy effects (see XI). 

X. PROPAGATION OF ELECTROMAGNETIC WAVES IN A TURBULENT 
CHANNEL FLOW 


An important application of the previous formalism is to the case of propagation of a 
laser beam through a turbulent channel flow medium. Of primary importance is, for 
example, the evaluation cf the degree of attenuation in the beam intensity I. Using the 
formalism developed by Hogge et al ., 33 the total beam intensity I, as a function of the 
coordinates x 0 and y 0 in the plane perpendicular to the direction of propagation, is given by 


0 f r r 

I(X °' yo) = (777 exp( " < 7 J ) J dx,dA 2 dy i dy 2 ex P[ C 0 ^) - g " >P/f x 0 (x r 
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where 


wfg = xf + xf + Vj + y 2 2 , 


(99) 


P 2 = (Xj-x 2 ) 2 +(y,-y 2 ) 2 , 


( 100 ) 


_ _ 2k 

P = ~- (101) 

Here. A is the wavelength, f the focal distance, and w, the spot size. In Eq. (98), (V(p) is 
the phase autocorrelation function and a 2 its value at p = 0. By separating the integrand 
in (98) so as to exhibit the coherent and incoherent parts of the total intensity, one ca" 
then evaluate the corresponding powers obtained by integrating the intensity over the 
variables x 0 and y 0 . A simple integration then gives the exact result 

? = P |nc /P t „,. 1 = 1 - eX P<- < ’«)' 002) 

which is often written as 


p - 1A 

1 + 

<t> 


(103) 


since in most 


cases is smaller than unity. The quantity aj = C^(0) is defined as 34 


9 9 


d 3 k 


L L 


COS 




cos 


MM-) 


COS 


u;(L-£") $ n (k.x\x") 


( 104 ) 
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where 


k 2 

= 2p ' K ~ = kj + 1<5 . k 2 = 4- a ; 2 (105) 

the (' and integrations are over the the photon optical path of length L. which is 

related to the full channel width 2D, x' = (2D/L)£', and x" = (2D/L)£". One of the 

integrations over optical path in (104) can be done by changing variables to sum and 

difference coordinates, (£'+£")/2 and ?/-£". We then assume that $ is a locallv 

n 

isotropic field, i.e., a function oniy of the sum coordinates, so that the integration over the 
difference coordinates can be done analytically. 

The function 4> n is defined as the spectral function of the square of the refractive 
index fluctuations, i.e., for mean flow in the Xj direction with the coordinate svstom 
oriented such that X 2 is the coordinate variable spanning the channel with 0 < x^ < A = 2D. 
we have 

<n ,2 (x^)> = J d 3 k $ n (k,x 2 ) , (106) 


where we have taken 
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n = 1 + n' = 1 + op, a = 7.9 * 10 ^g( A)/T . (107) 

Here the pressure p is measured in e.g.s. units and the temperature T in K. The slight 
dependence on the wavelength A is represented by the function g(A), where g(A=l//) = 1. 
g(A=0.5/<) = 1.02, and g(A=0.2/r) = 1.18. It then follows that 
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<p 2 (x 2 )> 


< n '*(x 2 )> = (ap) 2 = $ (k,x 2 ) d 3 k 

nr J n 


OOS) 


where p is the density in cgs units. Using (71), it follows that 


V k *x 2 ) = (opr n(k,x 2 ) . 


109) 


1 ) The quantity pa\j\ 


Using the equation of state and the sound speed c g 



8.32 * 10 7 

P 


P T, 



( 110 ) 


where k is Boltzmann's constant. N A is Avogadro's number, p the molecular weight 
(gm/mole), and 7 the ratio of specific heats, we derive using (107) 


poU 2 = 7 »M 2 , 


( 111 ) 


where the constant 7 * (p is in c.g.s. units) and the Mach number M are defined as 


7 * i 6.5696 7 jj, 



( 112 ) 


2 ) The results 

In Table V, we present the values of 10 2 P, Eq. ( 102 ), calculated using Eqs. (104). 
(109), and (94). The spectral function E(k) is calculated using both the C.ISS model and 
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the DIA. For each entry, ^and M, the corresponding U 0 , D (the channel half-width), 
and L (optical path) can be computed from the following relations (p in atmospheres and p 
in c.g.s. units) 


U 0 (cm/sec) = c g M = 10 3 M (tp Ip)* (113) 

D(cm) = , L(cm) = 2(1 + n 2 )* D(cm), (11-1) 

The values of the viscosity v (in cm 2 s _1 ), density p (in gm cm' 3 ), and the refractive index 
n r are given in Table V. A plot of 10 2 P versus the Mach number is displayed in Fig. (11) 
for A = 2, 3, and 5 mm. 

XI. ANISOTROPY * 


As mentioned at the end of Section IX, tve believe, on the basis of Kraichnan's 
work,” that the reason for the discrepancy between the calculated and measured values of 
the pressure fluctuations at the wall is the lack in our calculation of the inclusion of 
anisotropy effects induced by the mean flow. This is also believed to be the cause of the 
discrepancy between the calculated and observed one— dimensional energy spectra. 
Kraichnan presented a simple model for an elongated edoy structure by introducing a 
(constant) scale change in the direction of the mean flow. This simple model will be 
adopted here in order to judge what qualitative effect anisotropy would have on the 
one-dimensional energy spectrum. 

Following Kraichnan,” we introduce a transformation of scale defined by 
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Xj — ftXj , x 2 — , x 3 — X 3 , 


(1 loa) 


for the spatial coordinate variables where 1 is the streamwise direction and 2 is the 
direction across the channel, and 


I,/ _ 

K i - 


— k 2 , k 3 — k 3 , 


(115b) 


for the wavenumber variables, where a is a constant scale factor. The primed variables in 
(115), i.e., those for which the scales are stretched in the streamwise direction, are the 
physical variables. Since we only know E(k) for the isotropic case (in which there is no 
stretching), we derive a relationship between the velocity covariance tensors in the primed 
and unprimed systems. For homogeneous and isotropic turbulence, the velocity correlation 
function $jj(k) is given by (Ref. 15, Eq. 3.4.12) 




(116) 


where £jj is the Kronecker delta, E(k) is the energy spectrum, and 

<U J > = J d 3 k *„(k) . (117) 

The velocity covariance is a second rank covariant tensor. Applying the coordinate change 
(115b) to*ij(k) yields 


*n(k') = No^.fk) 


(118) 

* li (k') = No* lj (k), 


(119) 
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*ij(k') = N*jj(k), 


i,j / 1 , 


( 120 ) 

for the velocity covariance in the k' coordinate system, where tne normalization coefficient 
N has been added to these expressions to insure that the mean kinetic energy is identical to 
that of the original fiow. The new correlation tensor now describes a flow for which all 
lengths and velocities are elongated by a factor a in the x l direction. To calculate N. note 
that the kinetic energy (per unit mass) in the new coordinate system is 

) J *ii(k') d 3 k- = * ff^R,, + + R 33 ) , (121) 


where Rjj - / $ii(k) d k (no sum on the indices). To insure that the mean kinetic energy 
is identical to the original flow, (121) is equated to *(R n + + R 33 ). For isotropic R n . 

the resulting equation is solved for N and we find that 


a) One-Dimensional Spectra 


The one-dimensional energy spectrum in the physical coordinate system is given by 

E,(kj) = 2 JdV $> n (k') , (123) 

where k' 2 = k£ 2 + k 3 2 . After changing the variable of integration to k' 2 = a 2 k( 2 + h' 2 
and taking k'-k in the integral for ease of notation, the result for the one-dimensional 
energy spectrum in the physical variables is 
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( 124 ) 



Fig. (12) is a plot of the observed E,(lq) for R exp =30800 and that calculated using the 
DIA result. (Actually, we plot F^n), see Eq. (64).) The scale factor o is a free 
parameter and Fig. (12) displays the results for a = 1, 4, and 8. 

This procedure has provided a qualitative improvement over the a = 1 result. The 
turnover of the spectrum, originally occurring at n - 200 s -1 , has moved to frequencies 
(i.e., wavenumbers) beyond the lowest observed, in concert with the observations. 
Naturally, it is not to be expected that the procedure presented here would provide a 
perfect fit to the data since the mean flow may not affect all scale lengths in the same way. 
However, the results do suggest that accounting for anisotropy considerably improves the 
fitting of the one-dimensional spectra. 

b) Pressure Fluctuations 

For the anisotropic case, the calculation of the pressure fluctuations is similar to that 
done in Section IX. For the amplitude of the pressure fluctuations, we write 

a>i ^ )> = J d 3 k' n(k-,x;), (125) 


00 


n(k',x$) = 4sS^I Y «k$-l4 n) )$|^*»(k') fr;(k',x;) + Ifk-.x;) 


A'm4 

n=-oo 


( 126 ) 


r i\,' _ cosh k'(^A-x?) k' f, „xV 

1 l(k ’ X 2 } ~ sinh Ik'A + tjl 1 ~ 2 S“j 


4k 


sin kj x 2 jt" $ ^cos k^x^ (127) 
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r fir' V '\ - _ 2 _ f, oK' 2 lsinh K'fiA— Xi) . f. ,xV 
i(R ■ X;) ~EiSl 1 ~ 2 F^J cosh U- A l + E >- 2 r 


cos k^x 2 


4^ / 

4- j 77 2 ^sin k 2 x 2 , 


( 12S) 


with flx) the Dirac delta function, k^ n) = 2m/A, and k ' 2 = k ', 2 + k^ 2 . In (125)-(12S). 
the integration variables and position variable have been written as k' and x 2 , respectively, 
in order to emphasize that these are the physical coordinates, i.e., those stretched in the 
streamwise direction. Although we only know E(k) in terms of wavenumbers from which 
the stretching has been removed, we do know how to relate ^(k') to ^(k). Now change 
variables in the integral from k' to k. Working in the coordinate system defined by 


kj - k sin #cos 0 , k 2 = kcos#, k 3 = k sin # sin 0 , (129) 

the quantities k ' 2 and k ' 2 are given by 


*' 2 = k 2 sin 2 # (1 - cos 2 0 ) = k 2 f( 0 )sin 2 # , 

k /2 = k 2 (1 sin 2 # cos 2 0 ) = k 2 g(#, 0 ) . 

Using d 3 k' = d 3 k/o, we have from (125) 


^ ^ J d 3 k n(k,x 2 ;o) , 


(139) 


(131) 


(132) 


where 
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n=-oo 


E(k) sin 2 # cos 2 <t> k>5 

Ti[^[py 0 


X (r 2 (k',x 2 ) + r 2 (k',x 2 )) , (133) 

and T, and V 2 are given by (127) and (12. V Using (129) to write the delta function on k 2 
as a delta function on cos 0, the 9 integration in (132) can be done. Since the integrand is 
symmetric in the summation index n and the n = 0 term vanishes, we find for the 
amplitude of the pressure fluctuations 


<p 2 (xo;a)> _ 4N s3 

p 2 




00 


d<t> J dk 


E(k 


2i 7a 


COS 2 

“T[? 


[k A / 2 7rl 

T 

aLj 

n = l 


f2 Trnl 2 

1 [2 7T11 1 

2' 


' [ES"J 



g W~3) 


x [r-(k',x 2 ) + r^(k',x 2 )] , 


(131) 


where 


sin 2 0 n = 1 - 


2rn 

VET 


(135) 


Equation (134) was integrated numerically using the DIA spectrum for a = 4,8 and 
the results for Reynolds numbers of 12300, 30800, and 50000 are displayed in Figs. (13) and 
(14). Relative to the a = 1 case, the pressure at the wall is reduced by a factor of 
approximately 1/ a, i.e., with the pressure in units of /jU 2 , 

Pwall(a) « ^ Pwall(o^=l) • (136) 
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The pressures at mid— channel have undergone relatively small changes relative to the 
isotropic (a = 1) calculation, i.e., 


Proidcli(^) ~ Pmidchlu^l) • 


(137) 


It seems clear that the inclusion of the effects of anisotropy in the calculation of the 
amplitude of the pressure fluctuations has greatly improved our predictions of the wall 
pressures. 


c) Fraction of Scattered Power, P 


We have also applied the above procedure for the modeling of the effects of anisotropy 
on the attenuation of a laser beam propaga* : through a turbulent channel flow. With the 

•''.action of scattered power defined by (103) and (104), we integrated (104) over the 
physical variables, using the DIA spectrum and the expressions from above for the 
anisotropic pressure function. The scattered fraction P was calculated for channel widths 
A = 2, 3, and 5 mm at a variety of Mach numbers nd the results are presented in Table 
VI and displayed in Fig. (15) along with the experimental results 3 ' at these values of A./ 

|0 V. /-! rtr \ " * i ■ -s f ■ < 

The errors on the experimental results are roughly ±35%. The Mach number M u used in 

b 

Fig. (15) is the Mach number based on the bulk velocity, U . Using the first of (29) and 
the second of (30), we find that M b is related to the Mach number M (based on the 
midchannel velocity U 0 ) by 


M 


b 


M 


2.53 

9.268 Rexp 89 


(138) 


It was found that the model fit the data very well with a = 6 for A = 2 mm, a = 4 for A 
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= 3 mm, and a = 5 for A = 5 mm. These are the best fit results; however, within the 
errors, the data can be adequately fit for a = 5. It was found that over the range 1 < a < 
10, the scattered fraction P decreases with a roughly as 

Pta^P^lje"^ 1 )/ 3 . (139) 

XII. SUMMARY AND CONCLUSIONS 

By inspecting our results, it can be concluded that the GISS model and the DIA 
provide a reasonable first description of both turbulent bulk properties as well as properties 
that depend strongly on the large eddy part of the energy spectrum. Several reservations 
and directions for future work have to be discussed. 


First, the DIA yields larger values for the wall pressure and fraction of scattered 
power than the GISS model. This behavior can be understood by inspecting, for example. 
Eq. (94) which shows that the largest contribution to II(k), and ultimately to P, comes 
from die small k region of the turbulent energy spectrum. As one can see from Fig. 5, at 
low wave numbers the GISS E(k) is "skinnier" than the one derived from DIA. This means 
that the contribution to II(k) from the GISS model begins only at, say. kD =15, while in 
the DIA case the contribution begins earlier, i.e., at smaller wavenumbers. (The DIA 
energy spectrum extends in principle to zero wavenumber, but the integration over the 
angles in Eq. 92 and the presence of the delta function in Eq. 94 force the first contribution 
to begin at kD = ■k.) The physical reason behind the different low k behavior of the two 
spectral functions is known. The DIA model includes backscatter, i.e., energy transferred 
from the high k region into the low k region. By contrast, the present version of the GISS 
model is a cascade model with no backscatter as yet. While this limitation had been 
lecognized in earlier work, ’ it had never been evidenced as clearly as in the present 


50 


case, since our earlier works were primarily concerned with the description of turbulent 
bulk properties which are relatively insensitive to the small k behavior of the turbulent 
spectral function. Work to include backscatter in the GISS model is now in progress. 

A second comment concerns the role of shear and the way it has been treated in the 
present paper. The work of Tchen 11 ucu. demonstrated that the shear may play two roles, 
as a source of energy and as a force interacting with the eddies. While the first role is 
always present, the second becomes important only when the shear s is of the same order or 
larger than the eddy vorticity in a given wavenumber region. When this is the case, a 
resonance takes place which may be more important than the effects of non-linear transfer 
among the eddies, represented by the function T(k) in Eq. (2). When this type of 
resonance dominates over the other forces, Tchen's model predicts the existence of a k -1 
region in the energy spectrum which has indeed been observed. Our model does not include 
this resonance effect, the shear playing in fact only one role, that of a source of energy 
drained from the mean flow. This can clearly be seen by inspecting Eqs. (2)— (3) where the 
presence of shear is confined entirely to the left hand side of the equation in the rate n (k). 
whose form is given by Eq. (53). It is clearly seen that of the two terms, the first, Im(Q), 
is proportional to the shear, U'(y), while the other term represents a sink due to viseositv. 

A third comment refers to the choice of the growth rate. Orszag and Patera 35 
(hereafter referred to as O-P) have pointed out that the adoption of a Poiseuille flow with 
a central line velocity l T 0 (used to construct a Reynolds number which is then identified 
with the experimental value), leads to a growth rate that, being viscid in nature is 
naturally rather slow. In fact, the fastest instability occurs at around R=48000 and its 
maximum value (n g /n*) is only 0.0076. O-P discovered that there is a secondary, three 
dimensional instability that, catalyzed by the previous one, grows much faster than the 
original instability. The O-P discovery is an important one in many respects, but 
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particularly when one is concerned with the problem of the transition between laminarity 
and turbulence, since it explains in a natural way (i.e. without free parameters) 
experimentally important features that would otherwise be left unexplained, 36 the most 
important being the well known fact that linear analysis predicts a breakdown of 
laminarity at R exp =5772 while the experimental value is close to R ~2500. It is much 
less clear, however, whether the O-P secondary instability and its corresponding growth 
rate is directly relevant to our problem, where we have to deal with one important aspect 
that does not enter into the O-P problematics, namely the fact alluded to previously, that 
in the construction of n g for our problem we must somehow take into account the 
renormalizing effect of the presence of turbulence itself. The physical model that we have 
proposed in Section IV can be viewed as a way to renormalize the central line velocity due 
to the presence of turbulence. Another important aspect of the physical effect brought 
about by the existence of the O-P secondary waves is displayed in Figure 4 of the O-P 
paper, where it is shown that the growth rate versus Reynolds number curve reaches a 
maximum and then saturates rather than decreasing as it would in the case of a viscid 
instability. Two facts must be noticed. The maximum value of the O-P growth rate is 
around 0.1, which is the same as we have obtained using our renormalization (see Fig. 3). 
Secondly, as it is also clear from Fig. 3, in the region of Reynolds numbers of interest in 
this problem, the growth rates have similar maximum values, thus indicating that we are 
dealing with a region where the stabilizing effect of decreasing viscosity has not yet taken 
place. The numerical similarity between our results and those of O-P constitutes, 
however, no guarantee that the O-P mechanism (to generate a growth rate) and our 
method are physically equivalent. Since our results for the bulk properties are in general 
smaller than the experimental values, it is conceivable that the use of an O-P growth rate 
together with a renormalization of some of the physical parameters, could yield better 
results. Since, however, the solution of the O-P equation, which in turn requires the exact 
solution of the Navier-Stokes equations for the 2-dimensional instability, would be a 
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non— trivial addition to our problem, we decided in this first paper to adopt the physical 
argument discussed above. The use of the O-P growth rate mechanism is presently under 
investigation. 

In conclusion, while the use of theoretical models of turbulence has reproduced several 
properties of turbulent channel flow, limitations have also appeared. Work to include 
backscatter in the GISS model and to understand the role of the secondary instability is 
now in progress. 
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Table I 


R 

exp 

R r 

R m 

r i 

r 2 

10 4 

500 

8735 

12.55 

9.54 

2 * 10 4 

925 

17659 

21.00 

16.15 

3 * 10 4 

1331 

26633 

29.00 

22.20 


Table II 


The growth rate n g (k) vs. k from a solution of the Orr— Sommerfeld equation for four values 
of the Reynolds number. R (The unit n* is given by Eq. 52). 


R = 

exp 

12,300 

R - 

exp 

28,600 

R = 

exp 

30,800 

R = 

exp 

50,000 

kD 

10n g /n* 

kD 

10n g /n* 

kD 

10n s /n, 

kD 

10n g /n* 

8.631 

0.125 

10.533 

0.173 

10.700 

0.168 

12.74 

0.410 

9.081 

0.291 

11.265 

0.414 

11.256 

0.356 

14.51 

0.838 

9.588 

0.447 

12.113 

0.643 

11.886 

0.537 

16.57 

1.19 

10.144 

0.589 

13.053 

0.849 

12.578 

0.709 

18.85 

1.41 

10.739 

0.713 

14.068 

1.025 

13.323 

0.867 

21.33 

1.42 

11.369 

0.813 

15.147 

1.158 

14.112 

1.006 

24.0 

1.08 

12.031 

0.886 

16.285 

1.234 

14.940 

1.123 

27.1 

0.117 

12.721 

0.927 

17.479 

1.254 

15.804 

1.212 



13.441 

0.930 

18.732 

1.183 

16.702 

1.268 



14.189 

0.889 

20.053 

1.003 

7.634 

1.283 



14.972 

0.796 



18.599 

1.251 



15.794 

0.645 



19.603 

1.160 







20.650 

0.998 
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Table III 


ky 35 a function of k ± computed from Eq. (5S) for three values of the Reynolds number. 


R = 

exp 

k D 

i 

12.300 

k y D 

R = 

exp 

k D 

X 

30,800 

k y D 

R = 

exp 

k D 

X 

50,000 

k y D 

0.453 

8.62 

0.425 

13.3 

0.32 

12.7 

0.4S 

9.07 

0.485 

15.8 

0.36 

14.5 

0.507 

9.58 

0.505 

16.7 

0.40 

16.6 

0.533 

10.1 

0.525 

17.6 

0.44 

18.8 

0.56 

10.7 

0.545 

18.6 

0.48 

21.3 

0.587 

11.4 

0.565 

19.6 

0.52 

24.0 

0.613 

12.0 

0.585 

20.6 

0.56 

•» 

£ / . 1 

0.64 

12.7 

0.605 

21.7 





0.625 

22.9 
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Table IV 


R = 12.300 

exp 

GISS DIA EXP 


10 3 <U 2 > 0.920 1.14 

10 2 <U5>* 1.75 1.95 

10% 6.30 

10 5 f 4.70 5.43 

A x 0.192 0.092 

Ax 0.857 1.05 

10% 10.3 9.97 


3.6 


0.11 


R = 30,800 

exp 


R =50. (Hi!) 

exp 


GISS 

1.16 

1.96 

5.04 

7.12 

0.15 

0.60 

4.70 


DIA 

1.21 

2.01 

7.53 

0.051 

0.80 


EXP 

3.0 

3.62 

0.095 

0.63 


GISS 

1.25 

2.00 

4.30 

9.08 

0.141 

0.49 
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Table V 


Values of 10'P, Eq. (102), and full channel width A (in mm) 
for different Mach and Reynolds numbers. 


M 


R exp = 12300 


GISS 


DIA 


0.1 

5.71 

6.50 

X 

10- 4 

1.04 

x 10" 3 

0.114 

5.00 

8.49 

X 

10 -4 

1.36 

x 10~ 3 

0.190 

3.00 

2.36 

X 

10~ 3 

3.77 

x 10~ 3 

0.2 

2.86 

2.59 

X 

10“ 3 

4.16 

x 10 -3 

0.286 

2.00 

5.30 

X 

10 -3 

8.48 

x 10“ 3 

0.3 

1.90 

5.82 

X 

10~ 3 

9.35 

T 

o 

X 

0.4 

1.43 

1.04 

X 

10“ 2 

1.66 

X 10 3 

0.5 

1.14 

1.62 

X 

10 -2 

2.59 

X 10~ 2 

0.6 

0.952 

2.33 

X 

10' 2 

3.72 

C I 1 

o 

X 

0.7 

0.816 

3.23 

X 

10 -2 

5.07 

C I 1 

o 

rH 

X 
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= 30800 


M 


A 


GIS: > 


DIA 


0.1 

14.3 

9.04 x io -3 

1.81 x 10 

0.2 

7.15 

3.56 x 10 -3 

7.23 x 10 

0.286 

5.00 

7.38 x io -2 

0.148 

0.3 

4.77 

8.08 x io -2 

0.162 

0.4 

3.58 

0.144 

0.289 

0.477 

3.00 

0.205 

0.410 

0.5 

2.86 

0.225 

0.450 

0.6 

2.38 

0.324 

0.649 

0.7 

2.04 

0.438 

0.883 

0.715 

2.00 

0.459 

0.920 
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M 

A 

G1SS 

DIA 

. 

0.1 

23.2 

3.33 * 10 -2 

4.30 * 10 -2 


0.2 

11.6 

0.133 

0.183 

a 

0.3 

7.74 

0.297 

0.413 


0.4 

5.81 

0.530 

0.734 


0.465 

5.00 

0.714 

0.9S7 


0.5 

4.65 

0.827 

1.14 


0.6 

3.87 

1.19 

1.64 


0.7 

3.32 

1.59 

2.23 


0.774 

3.00 

1.97 

2.71 


1.16 

2.00 

4.35 

5.99 


Gas properties: 

C | > 

O 

X 

o 6 

if 

cm 2 sec -1 , p = 2.07* 10 -3 gm 

-3 
cm , 



p = 2 atm, 

p = 28 gm mole -1 . 



Optical properties: wavelength A = 

= 0.53 microns, refractive index 

n r = 1.53. 
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Table VI 


Values of 10 2 P for full channel width 
A = 2, 3, and 5 mm. 

A = 2 mm 


R exp 

M 

M b 


10 2 P 




a = 1 

co 

II 

<3 

II 

C5 

12,300 

0.286 

0.252 

8.48 « 10 -3 

2.80 x lCf 3 7.86 x io -4 

30,800 

0.715 

0.638 

0.920 

0.315 8.89 x io" 2 

50,000 

1.16 

1.04 

5.99 

2.25 0.662 


A = 3 mm 


R exp 

M 

M b 

10 2 P 

a=l a = 4 a=8 

12,300 

0.190 

0.168 

3.77 x 10" 3 1.24 x 10“ 3 3.49 x 10^ 

30,800 

0.477 

0.425 

0.410 0.140 3.95 x io~ 2 

50,000 

0.774 

0.694 

2.71 1.01 0.295 
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FIGURE CAPTIONS 


Figure 1. The expected general shape of the growth rate n.(k) vs. k. At largo 
wavenumbers, n (k) tends to -i/k 2 . 

Figure 2. The growth rate n g (k) vs. k in units of n*, Eq. (52), as solution of the 
Orr-Sommerfeld equation for F = 12300, 30800, and 50000. 

vAp 

Figure 3. The GISS turbulent energy spectral function F(k), in units of F*, Eq. (59), vs.k 
for R exp = 12300, 30800, and 50000. 

Figure 4. A comparison of the GISS and DIA spectral functions for R = 30800. (Same 

GXp 

units as in Fig. 4) 

Figure 5. The energy «(k) per unit mass and time, Eq. (7), in units of c*, Eq. (60), for 
several values of R exp . As one may note, the function saturates very rapidly, thus 
becoming independent of the wavenumber k. 

Figure 6. The turbulent viscosity i/Jk) vs. k, Eq. (8), in units of */*, Eq. (61), for several 
values of R 

exp 

Figure 7. Comparison of the theoretical one— dimensional spectral function, Eq. (63). 

vs. k t (full line), with the one measured by Hussain and Reynolds** for R = ^SGOO 

exp 

U 0 = 1350 cm sec 1 and D = 3.18 cm. 


66 


Figure 8. Comparison of the theoretical one-dimensional spectral function defined hv F.q. 

(64) with the one measured Laufer 1 , corresponding to R = 30800 U 0 =728 cm 

exp o 

sec and D = 6.35 cm. 

Figure 9. The value of the fluctuating pressure, in units of /?U 2 computed from Eq. (93). 
as a function of distance across the channel of width A=2D. The results arc- 
computed using the GISS model. 

Figure 10. Same as Figure 10, but computed using the DIA model. 

Figure 11. The quantity 10 2 P, Eq (102), versus M 4 for three values of the 
channel width A = 2, 3, 5 mm. 

Figure 12. Same as Fig. 8, but also for o= 4 and 8. 

Figure 13. Same as Fig. (10), but for a = 4. 

Figure 14. Same as Fig. (10), but for a = 8. 

Figure 15. The quantity 10 2 P versus the fourth power of the bulk Mach number, Mj} (see 

text), for three values of channel width A = 2, 3. and 5 mm. The best fit value of a is 

shown. The experimental data are represented as 7 for the 5 mm results, * for the 3 
mm results, and o for the 2 mm results. 


67 













Fig. 5 




Fig. 6 



E,(kj) 



Fig. 7 












F,2(n)sec 



Fig. 12 








Fig. 15 






